Quantum information processing using strongly-dipolar coupled nuclear spins 
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Dipolar coupled homonuclear spins present challenging, yet useful systems for quantum informa- 
tion processing. In such systems, eigenbasis of the system Hamiltonian is the appropriate com- 
putational basis and coherent control can be achieved by specially designed strongly modulating 
pulses. In this letter we describe the first experimental implementation of the quantum algorithm 
for numerical gradient estimation on the eigenbasis of a four spin system. 
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An important issue in experimental quantum informa- 
tion processing (QIP) is achieving coherent control while 
increasing the number of qubits 0, 0- In all existing 
implementations of quantum computation, the execution 
time is limited by durations of nonlocal gates. In the case 
of nuclear magnetic resonance (NMR) implementations 
of QIP, where qubits are formed by mutually coupled 
spin 1/2 nuclei, most of the existing implementations 
used isotropic liquids. In these systems, the durations 
of nonlocal gates are limited by the strength of (indi- 
rect) scalar couplings, which typically are of the order of 
10-100 Hz. 

Apart from the scalar couplings, nuclear spins also in- 
teract through magnetic dipole-dipole couplings, which 
are 2-3 orders of magnitude stronger than the scalar cou- 
plings, and would therefore allow significantly faster gate 
operations. In the case of isotropic liquids, the rapid 
molecular reorientation averages the dipolar interactions 
to zero. On the other hand, in oriented systems, the 
dipolar interactions survive and are therefore potentially 
better candidates for NMR-QIP . 

The Hamiltonian for a dipolar coupled n-spin system 

is, 

n n 

H = J2^4+ E 27^(31^ -F -I fc ) (1) 
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where ujj are the chemical shifts, Djk are the dipolar 
coupling constants and I J Z are z-components of the spin 
angular momentum operators V [^|. The first term 
corresponds to the Zeeman interaction and the second 
term describes the dipolar interaction. When \Djk\ <C 
\u)j— Wk\, one usually invokes the 'weak-coupling' approx- 
imation 0, [j| , where one ignores the off-diagonal parts 
of the Hamiltonian. Then the Zeeman product states 
are the eigenstates of the full Hamiltonian and individ- 
ual spins can be conveniently treated as qubits. Most of 
the experiments in NMR-QIP have so far been carried 
out on such systems [2|. 

However to maximize the execution speed of our quan- 
tum processor it is desirable to use stronger couplings, 
including |D^| > \uij — u>k\- In this case, we have to use 
the complete Hamiltonian (1), where the Zeeman and the 
coupling parts do not commute and not all eigenstates 
are Zeeman product states but linear combinations of 



them. Unlike in the case of a weakly coupled spin sys- 
tem, individual spins of a strongly coupled system loose 
addressability. In such cases the eigenstates of the full 
Hamiltonian can form individually controllable subsys- 
tems and then the eigenbasis becomes the natural and 
accurate choice for the computational basis. 

Though strongly dipolar coupled spin-systems have 
been suggested for QIP earlier , so far it had not been 
possible to implement general unitary gate operations 
that can be applied to arbitrary initial conditions. Here 
we show how quantum computation can be realized on 
the eigenbasis by efficient coherent control techniques. As 
a specific system, we use a strongly dipolar coupled four- 
spin system partially oriented in a liquid crystal. Such 
systems have certain key merits. Unlike in the liquid 
state systems, the intramolecular dipolar couplings are 
not averaged out completely, but are only scaled down 
by the order parameter of the solute molecules oriented 
in the liquid crystal 0, @ ■ Since intermolecular interac- 
tions are averaged out (unlike the crystalline systems), 
liquid crystalline systems provide well defined quantum 
registers with low decoherence rates. We achieve high- 
fidelity coherent control with the help of specially de- 
signed strongly modulating pulses (SMPs). 

To demonstrate these techniques, we use an interesting 
algorithm suggested by S. P. Jordan |9(: the quantum al- 
gorithm for numerical gradient estimation (QNGE). Be- 
fore we discuss our implementation, we summarize the 
theoretical description of QNGE 0. The gradient of an 
one-dimensional real function / over a small real range 
I is written as V/ = {f(l/2) - f(-l/2)}/l. Thus classi- 
cally two function evaluations are necessary to estimate 
the gradient of a one-dimensional function and a mini- 
mum of d + 1 function evaluations are necessary for a 
d-dimensional function. On the other hand, QNGE re- 
quires only one function evaluation independent of the 
dimension of the function. Here the function / is en- 
coded in an n-qubit input register. An ancilla register 
is also required whose size (no qubits) depends on the 
maximum possible value of the gradient. 

In QIP, numbers are represented in binary form. To 
encode the real number x in the input register, we have to 
convert it into a nonnegative integer <5 £ {0, 1, • ■ • , TV— 1}, 
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FIG. 1: Circuit diagram for QNGE. 

where N = 2™. The encoding is defined by 
I f, N-l" 



N - 1 



(2) 



The circuit diagram for the quantum algorithm is shown 
in Figure 2] Initially the ancilla register (an) is set to 1 
and the input register (in) is set to 0. The inverse quan- 
tum Fourier transform (IQFT) prepares a plane wave 
state on the ancilla register and and the Hadamard trans- 
form (H) prepares a uniform superposition on the input 
register [lOj: 
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The ancilla register is now in an eigenstate of the ad- 
dition modulo iVo- On applying oracle Uf : \an,in) — > 
\sf(x) ®n an, in), where s is a scaling factor, the total 
state becomes 
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For a small x, f(x) « f(0) + xVf. Substituting for x 
from equation |2 and ignoring the global phase, the input 
register reduces to, 



g JVq AT - 1 v J 
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The scaling factor s is set to maximize the precision i.e., 

s = N (N -l)/Nl. (3) 

Now a phase estimation is carried out with the quantum 
Fourier transform (QFT) on the input register 
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Measuring now the input registers in the computational 
basis gives an estimation of V/. 
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FIG. 2: The SMP for QNGE: (top) amplitude vs time, (mid- 
dle) phase vs time, and (bottom) fidelity vs RF inhomogeneity 
(RFI, 0.9 to 1.05 of the ideal field) and static field inhomo- 
geneity (SFI, -5 Hz to 5 Hz). The 30-segment SMP has an 
average fidelity F' avg = 0.981. 



In our 4-qubit case, we select two qubits as ancillas 
and the other two as input qubits i.e., no = n = 2, 
N = N = 4, and S € {|00), |01), 1 10) , 1 11)}. From 
equation x G {-1/2,-1/6,1/6,1/2} for I = 1. Let 
us consider an example function f(x) £ {0,2/3,4/3,2}, 
which has a gradient V/ = 2. Using equation [3] we ob- 
tain the scaling factor s = 3 so that sf(x) £ {0,2,4,6}. 
For S £ {|00),|10)}, s/(x)eAr remains identity. For 
S £ {|01), |11)}, sf(x)<$N adds 2 to the ancilla, i.e., flips 
the first ancilla qubit. Therefore, for this example, the 
oracle is a CNOT(ani,wi2) gate i.e., a NOT operation 
on the first ancilla qubit controlled by the second input 
qubit. The propagator for the entire algorithm is there- 
fore, 
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We design a single SMP corresponding to this entire op- 
eration on the input states. 

An SMP is a cascade of radio-frequency (RF) pulses 
numerically calculated based on the precise knowledge 
of the internal Hamiltonian of the qubit system and the 
target propagator 0,0]. Given a target operator Ut, 
Fortunato et al 0] have described numerically searching 
an SMP propagator C/smp based on four RF parameters 
for each segment k: duration (t^), amplitude (w^), 

phase (cf>^) and frequency (uip). We found that it is 
useful to add one more degree of freedom: after each 

(k) 

pulse segment, we introduce a variable delay r d . The 
delays are computationally inexpensive to optimize, easy 
and accurate to implement and make designing non-local 
gates easier. All the parameters are determined so as to 
maximize the fidelity 

F = trace [U^ 1 ■ C/ S mp] /M 
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FIG. 3: Molecular structure of l-Chloro-2-iodobenzene (top- 
left), the elements of diagonalized Hamiltonian (top- right), 
and the energy level diagram with prominent transitions (bot- 
tom). The eigenstates are labeled as shown. 



where J7smp is the propagator for the SMP, M is the di- 
mension of the operators and s is the number of segments. 
A MATLAB package has been developed which uses the 
Nelder-Mead simplex algorithm as the maximization rou- 
tine 01 . The search constraints can be simplified if the 
input state is definitely known. The fidelity of an SMP 
specific to a known initial state p ln can be written as, 



F' = 



trace [p T ■ p SM p] 
>/trace [p%] ■ trace [p 2 SMP ] ' 



(4) 
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where p T = U T ■ p m ■ U^ 1 and p S MP = U S mp ■ Ph 
The SMPs are made robust against the spatial inho- 
mogeneous distributions of RF amplitudes and of static 
fields (Aa>j) by maximizing an average fidelity, F' avg — 
Y li F'(uj l A,Auj l ) FigureEJshows a single SMP per- 
forming the full QNGE on the specific four-qubit system. 
Though it is possible to decompose the target propagator 
into several SMPs each corresponding to one or two-qubit 
gates, it is more efficient, at least for small spin systems, 
to design and execute a single robust SMP implementing 
the entire algorithm. 

As the quantum register, we used the four X H spins 
of l-Chloro-2-iodobenzene (CIB; Figure EJ) (purchased 
from Sigma Aldrich®) oriented in liquid crystal ZLF 
1132 (purchased from Merck®) forming a lOmM solu- 
tion. All experiments were carried out on a 500 MHz 
Bruker Avance spectrometer at 300 K. 
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FIG. 4: 1 H spectra of CIB obtained by applying a small 
angle pulse (9 = 3°) to the equilibrium state (bottom), 
|0100) (0100| — 10000) (0000 POPS (middle), and QNGE (top). 
The transitions numbers correspond to those in Figure^] The 
pulse sequences for POPS and QNGE are shown in insets, 
(by applying a transition-selective Gaussian n pulse of 10 ms 
duration) Pulsed field gradients (G z ) are used to destroy co- 
herence after POPS and after the SMP. In the QNGE exper- 
iment, the residual zero-quantum coherence are de-phased by 
a 32-step average over a random delay r (between and 10 



Figure 0] shows the 1 H NMR spectrum of the partially 
oriented CIB. The line widths of the various transitions 
range from 1.7 Hz to 4.0 Hz indicating coherence times 
(T 2 * relaxation times) between 1.8 s and 0.7 s. The co- 
herence times are sufficiently long to ignore relaxation 
effects in the design of the SMP. 

The procedure for analyzing the NMR spectra of par- 
tially oriented systems has been well studied yj 13 ■ We 
have developed a numerical procedure to iteratively de- 
termine the system Hamiltonian from its spectrum and 
a guess Hamiltonian ^5|. The 37 strongest transitions 
of the CIB spectrum were used and a unique fit was ob- 
tained. The mean frequency and intensity errors between 
the experimental and the calculated spectra are less than 
0.1 Hz and 6% respectively. The elements of the diag- 
onalized system Hamiltonian and the corresponding en- 
ergy level diagram are shown in Figure 

In NMR-QIP, the initial states are not pure states but 
are pseudopure states that are isomorphic to pure states. 
The pseudopure states differ from the pure states by a 
uniform background population on all states. It is easier, 
however, to prepare a pair of pseudopure states (POPS) 
H3- We prepared the pair |0100)(0100| - |0000)(0000|. 
The first term represents the desired initial state; the ad- 
ditional second part does not interfere with the QNGE 
experiment, because the operation IQFT, when acted 
on 1 00), creates a uniform superposition of the ancilla 
qubits. Such a state is invariant under the f(x)®N °P~ 
eration and therefore the output state corresponding to 
1 0000) (0000 1 is independent of the oracle U f . The inset 
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FIG. 5: Bar-plots showing the diagonal elements of the den- 
sity matrix in the eigenbasis: equilibrium (top left), simulated 
POPS (middle left), experimental POPS (bottom left), theo- 
retical output of QNGE (top right), simulated output (with 
the SMP shown in Figure middle right), and experimental 
output of the quantum algorithm (bottom right). 



at the center of Figure 0] shows the pulse sequence for 
preparing POPS and the corresponding spectrum. The 
POPS spectrum is obtained by linear detection pulse af- 
ter inverting transition 2 and subtracting from the re- 
sulting spectrum the linearly detected spectrum of the 
equilibrium state. 

Since we use eigenbasis as the computational basis, the 
projective measurement on to the computational basis is 
equivalent to measurement of populations after destroy- 
ing the coherence. First we dephase the coherence by 
using a pulsed field gradient (PFG). A PFG does not ef- 
ficiently dephase homonuclear zero-quantum coherence. 
Therefore we use a random delay (in between and 10 
ms) after the PFG and average over several (32) tran- 
sients. The populations are then measured using a linear 
detection (small flip-angle) pulse 0. Using the normal- 
ization condition, there are 15 unknowns for the 16 eigen- 
states. The results of the diagonal-tomography obtained 
by the mean of three sets each of 15 linearly independent 
transitions are shown in Figure [SJ 

After the projective measurement in the eigenbasis at 
the end of the quantum algorithm, the input qubits are 
in the state 1 10) encoding the gradient V/ = 2, while 



the ancilla qubits have equal probability in all possible 
states. Therefore, the theoretical output diagonal state of 
the combined system is I an <gi (|10)(10| - |00)(00|), where 
I an is the identity operator for the ancilla qubits and the 
two parts in the parenthesis correspond to the two parts 
of the POPS. 

The diagonal correlation C between the theoretical 
density matrix (px) and the experimental density matrix 
(pe), is defined as 



C 



trace[|p T | ■ |pe| 



\/trace[|p T | 2 ] • trace[|p E | 



(5) 



where || denotes extraction of the diagonal part. The 
average diagonal correlations were 0.999 and 0.979 for 
the POPS and the result of the QNGE, respectively. The 
lower value of the latter may be attributed to decoherence 
and spectrometer nonlincarities. 

In conclusion, we have demonstrated quantum compu- 
tation on the eigenbasis of system Hamiltonian using co- 
herent control techniques. As an example, we described 
the first implementation of Jordan's algorithm for nu- 
merical gradient estimation. Compared to the usual ap- 
proach using weakly coupled systems, the present method 
using strongly dipolar coupled systems yields signifi- 
cantly faster execution times and is therefore less suscep- 
tible to decoherence effects. From molecular and spec- 
troscopic considerations a combination of homo and het- 
cronuclcar spins in cither liquid crystalline or molecular 
single crystalline environments is a natural way to build 
larger qubit-systems (with > 10 qubits). The coherent 
control of strongly dipolar coupled systems then becomes 
important and the present work is the first step in this di- 
rection. We believe that the coherent control techniques 
demonstrated here for dipolar coupled nuclear spins will 
turn out to be essential also for other solid-state imple- 
mentations of quantum information processing. 
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